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Abstract. We introduce a family of exactly-solvable two-dimensional Hamiltonians 
whose wave functions are given in terms of Laguerre and exceptional Jacobi 
polynomials. The Hamiltonians contain purely quantum terms which vanish in the 
classical limit leaving only a previously known family of superintegrable systems. 
Additional, higher-order integrals of motion are constructed from ladder operators 
for the considered orthogonal polynomials proving the quantum system to be 
CN ! superintegrable. 
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1. Introduction 

• 

The connection between classical families of orthogonal polynomials and exactly- 
solvable, integrable and superintegrable systems is well-known. In this paper, we show 
that the same connection applies to the recently discovered families of exceptional 
orthogonal polynomials In particular, we demonstrate the existence of an infinite 
family of Hamiltonians which are both superintegrable and exactly-solvable and whose 
wavefunctions are composed of the product of classical Laguerre and exceptional Jacobi 
polynomials. Furthermore, in the classical limit, the Hamiltonian reduces to the 
celebrated Tremblay-Turbiner- Winternitz system [2j |3] . 

A system is said to be superintegrable if it admits more integrals of motion than 
degrees of freedom. For the purpose of this article, we consider a Hamiltonian on the 
two-dimensional Euclidean plane 

H=^(pl+pl) + V(x u x 2 ) (1.1) 
which admits two addition integrals of motion 

L a = ^2 fa,jk(xi,x 2 )p{p2, a = 1, 2, (1.2) 

0<j+k<n 
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where the Pi are the components of the momenta conjugate to Xj and are taken to be, 
in the quantum system, 

d 

Pi = ~^~q^: C 1 - 3 ) 

As indicated in (jl.2p . the additional integrals will be assumed to be polynomial in 
the momenta and the degree of the system is said to be n, the highest order of 
the integrals. In classical mechanics the three functions L 2 , H} are assumed 
to be functionally independent and in quantum mechanics they are assumed to be 
algebraically independent. While the integrals, L a , Lie or Poisson commute with the 
Hamiltonian H, they do not commute with each other and so generate an algebra which 
usually closes to form a polynomial algebra [U El El [7] . 

The study of superintegrable systems began with second-order superintegrability 
[H El HOI HU EG] , including the best-known examples of the harmonic oscillator [13l [14] 
and the Kepler-Coulomb system [151 HE]- Recently, new families of superintegrable 
systems have been discovered with integrals of arbitrary order. The first was discovered 
by Tremblay, Turbiner and Winternitz [2j E] 

*™ + ^ + k ; (-^ + . (1.4, 

The discovery of these new families has lead to much new research on the discovery and 
treatment of such superintegrable systems, see e.g. [T71 [HI EEl 120] . 

The connection between superintegrable systems and orthogonal polynomials is 
most obviously evident in the conjecture that all superintegrable systems are exactly- 
solvable [2T]. Recall, a system is said to be exactly solvable if its energy values 
can be calculated algebraically and the wave functions can be written in terms of 
orthogonal polynomials multiplied by the ground state [221 E31 El]. This connection 
was recently exploited by Kalnins, Kress and Miller who made use of ladder operators 
for orthogonal polynomials to construct additional integrals of motion, thus proving 
the superintegrability of several families of superintegrable systems, including the 
TTW system [25]. This method has been successfully been applied to families of 
superintegrable systems with reflections in the potential |26j as well as scalar potentials 
defined on pseudo- Euclidean space [27J. Ladder operators associated with shape 
invariance have also been used to construct infinite families of superintegrable systems 

[2111211 ED]. 

The purpose of this article is to extend this analysis to exactly-solvable systems 
whose wave functions are expressible in terms of exception orthogonal polynomials. 
Exceptional orthogonal polynomials are eigenf unctions of Sturm-Liouville equations 
which generalize the classical families of orthogonal polynomials in the sense that, unlike 
the case of classical polynomials, the families of exceptional orthogonal polynomials 
admit gaps in their degree [H [31]. The Sturm-Liouville equations associated with 
such polynomials have been recently the subject of intense research including their 
connection with shape invariance, see e.g. [321 [331 EH ESI EHl E7J, and their connection 
with Darboux-Crumb transformations, see e.g. [381 EHl HQ]- 
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The plan of the paper is as follows. In section 2, a new exactly-solvable two- 
dimensional Hamiltonian is introduced and its wavefunctions and spectrum are found 
as well as its classical limit. In section 3 the Hamiltonian is proven to be superintegrable 
for rational k by direct construction of the integrals of motion. Section 4 gives the k = 1 
example and explores the connection with its classical limit. Section 5 is comprised of 
concluding remarks. 



2. An exactly-solvable Hamiltonian 

Consider a generalization of the TTW Hamiltonian (11.41) given in polar coordinates by 

H - *A I l ^ I fc V ft2 -i | | 4(l + frcos(2fc0)) \ 

* 2 A + 2 U + 2r* Un 2 (fc0) cosW) + (b + cos(2fc0)) 2 J 1 ' 

where 

_ (3 + a 
(3 — a 

The Schrodinger equation associated with this Hamiltonian ( 12. ip 

H k ^ — E^/ = 0. (2.2) 
separates in polar coordinates as 

$ = $(<j>)R(r) (2.3) 



with 



~d r rd r + X -ur 2 + ^--E) R(r) = (2.4) 



V, q2 -I i | 4(l + &cos(2fc0)) 2) 

P^ + sin 2 (^) + cosW) + (6 + cosW))2 W)- 0.(2.5) 

Note that the radial equation is exactly that of a two-dimensional oscillator and the 
angular part is a deformation of a Darboux-Poschl- Teller potential |41[ |4"2]. 

In fact, the angular part of the Hamiltonian (12.11) is the shape- invariant Hamiltonian 
introduced in [38] whose eigenfunctions can be written in terms of the X\ Jacobi 
polynomials. The change of variables 

$(0) = X n (x), x = cos(2A;0), n > 1, (2.6) 

transforms (12. 5p into 

(GxT^G- 1 - A 2 ) X n (x) = (2.7) 

where 

(x — b) 

T a,p = 4(x 2 _ l)d 2 + 4(/3- a)(l-fa;) + _ + + + 2 ^ 

— X 
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The operator T a ^ is equivalent to the eigenvalue operator for the X\ exceptional 
Jacobi polynomials [lj. The eigenvectors of (I2.9p are given by the exceptional Jacobi 

- — -a,/8 

polynomials P n , 

= -\{x - b)P»Z(x) + 2n _ 2 \ a + p [bPSAix) - P^ 2 (x)] , n > 1, (2.10) 

where P^ ,l3 (x) are the standard Jacobi polynomials, [38]. The eigenvalue equation (12.71) 
has eigenvalue 

A 2 = A 2 n = (2n - 1 + a + (3) 2 , n > 1. (2.11) 
For (12. 4p , the change of variables 

R(r) = Y^(y), y = ur 2 (2.12) 
and conjugation by the ground state 

Gy=y A n /2 e -y/2^ ( 213) 

transforms (\2A\i to 

yd 2 Y + (1 + /cA n - y)d y Y + ^Y = 0. (2.14) 

y 4a; 

The solutions of (12. 14f) are given in terms of Laguerre polynomials 

Y^ = G y Lt{y), (2.15) 

whenever the energy is quantized as 

E m>n = w(2m+ k(2n + a + /? - 1) + 1) . (2.16) 

Unless needed to avoid confusion, the subscripts on A and E will be dropped in what 
follows. 

2.1. The classical limit 

It is interesting to observe that in the classical limit, this system reduces to the TTW 
system (2]. In fact, we re-introduce the parameter h by multiplying the system by a 
factor of h 2 . The corresponding Schrodinger equation is equivalent to the h~ = 1 case, 
with an appropriate scaling of the energy. The Hamiltonian becomes 

H- ft^A | V.2 r 2 , ^V « 2 -i i P 2 -l i 4(l + fecos(2fc0)) \ 
^-" n 2 A+ 2^ r + (6 + cosW)) 2 J' (2 ' 17) 

We renormalized the parameters so that the potential is not annihilated in the classical 
limit {H — > 0) by taking 

a = ha, (3 = h[3, Q = hw. (2.18) 

The new Hamiltonian becomes 



1 1 P 

fe 2 2 2r 2 



-2_/^ /3 2 -^ 4 ^ 2 ( x + bcos(2k(f)) 
sin 2 (A;0) + cos 2 (^) H ^+^(2^" 



(2.19) 
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where 

£ = £+a lim6 = 6. (2.20) 

p-a h ~>° 

In this form, it is easy to see that the classical limit of the system (12. 19[) is exactly the 
TTW system and respects the requirement for bounded trajectories. Namely, in this 
notation, all bounded trajectories satisfy Q 2 > 0, a 2 > 0, f3 2 > 0. 

3. Superintegrability 

It is immediately obvious that the Hamiltonian (12. ip is integrable with first integral 
associated with separation of variables and given by 

1 d 2 ( a 2 -\ (3 2 -\ 4(l + fecos(2A;</>)) \ 
1 " k 2 d9 2 Vsin 2 (fc0) + cosW) (b + cos(2A;0)) 2 J' 1 ' ' 

In this section, we will show that the Hamiltonian is also superintegrable by constructing 
two additional integrals of motion for the Hamiltonian (12. ip using ladder operators for 
the orthogonal polynomials of the wavef unctions, as in [23] . 

The key to the method is to utilize ladder operators which transform the wave 
functions but leave the energy fixed, i.e. automorphisms on the energy eigenspaces. To 
this end, assume k = p/q, then the transformations 

m — )■ m + p, n — > n — q (3-2) 

m — >m — p, n^-n + q, (3.3) 

do not change the energy (12.161) . 

3.1. Ladder operators for the exceptional Jacobi polynomials 

Ladder operators for the exceptional Jacobi polynomials can be constructed from ladder 
operators for the Jacobi polynomials (see e.g.[43j) 

(1 -x 2 )(2n + a + (3) „ n ((2n + a + (3)x + a - (3 + 2) 
£n = d x (3.4) 

(1 -x)(2n + a + (3 + 2) ^ | (n + a + (3 + 1) ((2n + a + 0)x + a-p + 2) 
K-n ~ 2 x 2 

and the "forward" and "backward" operators for the exceptional Jacobi polynomials 

m 

J r =(*- 1 >(*+Hf)4 + <«- 1 )'(*+ a ^) w 
B °- .^°t^ " Ul)Ufl ' (3 ' 7) 

whose actions are 

FP^P-\ X ) = {2n-2 + 2a)P:f 1 (x) 



Exceptional Superintegrability 6 

The corresponding ladder operators for the exceptional Jacobi polynomials are 
defined as 

Cx, n = FoC n oB, n hn = FoTZ n oB. (3.8) 

Conjugating by the ground state, G x as in (12. 8p . gives ladder operator for the angular 
component of the wave function 

J-,n = G x L\^ n G x , J+,n = G x lZ\^ n G x l , (3-9) 

J_„X n (aO = -(n + a)(n + a - 2)(n + /3)(n + /3 - 2)X n _ 1 (x) (3.10) 

J+, n X n (x) = -(n)(n + /3)(n + a)(n + a + /3)X n+1 (x). (3.11) 
The repeated application of the operators J± jn are defined as 

J±,n = J±,n±(q-1) ° • • • ° J±,n±l ° ^±,n, (3.12) 

with action on the basis as 



J q _ >n X n = - a) q {-n -a + 2) q {-n - f3) q (-n -(3 + 2) q X n _ q 

Jl >n X n = (-l) q (n) q (n + 0) q {n + a) q (n + a + (3) q X n+q . 



(3.13) 



3.2. Ladder operators for the Laguerre polynomials 

Based on the ladder operators for the Laguerre polynomials, again see e.g. [43J, ladder 
operators for the functions Y^ A (y) can be constructed as 

T F kA 1 

K± kA , E Y* A (y) = (l±kA)d y - — T ^j(l + kA) Y* A (y) (3.14) 

= k ± Y^f (y), 

where 

fc+ = -l, k- = -(m + l)(m + kA), (3.15) 

and E and A take the quantized values as above (I2.16p . (12. lip 

E = 2u[2m + kA + 1], A = 2n + a + f3-l 

To obtain the desired shift in the quantum numbers (13.21) . the p-fold composition of K 
is defined with the corresponding value of kA shifted in each successive application 

KkA,E = KkA+2(p-l),E ■ ■ ■ KkA+2,EKkA,E (3.16) 
K-kA,E = K-(kA-2(p-l)),E ■ ■ ■ K-{kA-2),EK-kA,E- (3-17) 

The action of these operators on the wave functions is given by 

Ty-p ykA ( _i \pykA+2p 

ly kA,E 1 m \ L J 1 m—p fo io'i 

K P . kA , E Yi A = (-l)'(m + l) P {kA + m- P + l) p Y* A f p . { ' ] 

It is important to note that although the quantity E is a function of m and n, it is 
unchanged by the operation m — > mil, and kA — > kA =p 2. Hence, the energy E 
remains fixed in each successive applications of the operator. 
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3.3. Quantum-number independent integrals of motion 

Combing the two sets of ladder operators in subsections 13.11 and 13.21 operators can 
be constructed which transform within fixed energy eigenspaces. The corresponding 
operators are 

n + = ^kA,E^+,n^ n - = ^-kA,E^-,ni (3.19) 

which depend on the quantum numbers m and n and fix the energy E, so in fact 
[H ± ,#]* m , n = 0. 

In order to obtain differential operators which commute with the Hamiltonian for all 
values of m and n, these quantum numbers must be removed from the operators. To 
do this, first, the energy E is removed by moving the constant E to the right and 
the replacing it with H. See the derivation in the appendix of [26J for an explicit 
representation of these operators. To remove the quantum number n, note that under 
the transformation n — > —n — a — j5 + 1, the operators C n and lZ n (I3.4p .( l3~5l) transform 
as 

C n — > 7Z n , £ n +e — > 7tn-i- (3.20) 

This action transfers to the operators £i jn and lZ\ )n . Similarly, the action n — > 
—n — a — + 1 sends A — > —A and so interchanges the raising and lowering operators 
for the Laguerre polynomials, after the energy E has been replaced by H. Thus, the 
operators 

H 2 = S + + H_, ~ 3 = 1 - (5 + - 5_) , (3.21) 

2n + a + p — 1 

are polynomial in A 2 and hence the following operators 

L2 = ((H 2 ) E=H ) A 2 = ±i , L 3 = ((S 3 ) E=H ) A% u (3.22) 

fc 2 

are independent of the quantum numbers m and n and commute with the Hamiltonian 
H. Thus, L 2 and L 3 commute with the Hamiltonian on the basis and hence commute as 
operators based on a combination of writing the commutator, [Li, H], as a second-order 
operator with coefficients that depend on H and L\ [20] and using a standard Wronskian 
argument about the separated solutions [25] . 

The algebra relations of the integrals can be obtained directly from the expansion 
coefficients ( 13.131) . ( 13. 18}) . including the fact that neither L 2 nor L3 commutes with L\ 
so the operators are algebraically independent. 



4. Conclusions 



In this paper, we have introduced an exactly-solvable system whose wavefunctions 
are given in terms of a product of Laguerre and exceptional Jacobi polynomials. 
By construction, the Hamiltonian (12. ip is integrable and, in addition, it admits 
higher, independent integrals of motion making it superintegrable. These higher-order 
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integrals are constructed directly from ladder operators for the considered orthogonal 
polynomials. In particular, the ladder operators for the exceptional Jacobi polynomials 
are constructed from the "forward" and "backward" operators (13.61) . (13. 7h composed 
with the ladder operators for the classical Jacobi polynomials (13 .4p . (13. 5p . 

It should be noted than this method can be directly applied to prove the 
superintegrability of the following Hamiltonian, with a Kepler- Coulomb type potential 
in the radial variable, 



In fact, this Hamiltonian is related to the Hamiltonian (12.1 ft via coupling constant 
metamorphosis, see [T8 | l4"4" t H5]. 

This method of constructing Hamiltonians and their integrals of motion can be 
extended in a straightforward manner to other families of exceptional polynomials. Most 
immediately, the angular part of the Hamiltonian given above (12. II) can be replaced 
by any of the infinite families of one-dimensional Hamiltonians for the X% exceptional 
polynomials [331 EH]. Additionally, other families of Hamiltonians, say separable in 
Cartesian coordinates, can be obtained in a similar way from the Sturm-Liouville 
equations for other exceptional polynomials, e.g. extensions of the singular harmonic 
oscillator via exceptional Laguerre polynomials. These systems will be treated in future 
work. 
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